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ABSTRACT 

Dissipationless (gas-free or "dry") mergers have been suggested to play a major role in the formation 
and evolution of early-type galaxies, particularly in growing their mass and size without altering their 
stellar populations. We perform a new test of the dry merger hypothesis by comparing N-body 
simulations of realistic systems to empirical constraints provided by recent studies of lens early-type 
galaxies. We find that major and minor dry mergers: i) preserve the nearly isothermal structure 
(ptot r~ 2 ) of early-type galaxies within the observed scatter; ii) do not change more than the 
observed scatter the ratio between total mass M and "virial" mass R c a 2 2 /^G (where R c is the half-light 
radius and cr C 2 the projected velocity dispersion); iii) increase strongly galaxy sizes (R c cx M °- 85±0 - 17 ) 
and weakly velocity dispersions (<7 e 2 oc JVf 06 ± - 08 ) with mass, thus moving galaxies away from the 
local observed M-R e and M-a e 2 relations; iv) introduce substantial scatter in the M-R e and M- 
a e 2 relations. Our findings imply that — unless there is a high degree of fine tuning of the mix of 
progenitors and types of interactions — present-day massive early-type galaxies cannot have assembled 
more than ~ 50% of their mass, and increased their size by more than a factor ~ 1.8, via dry merging. 
Subject headings: galaxies: elliptical and lenticular, cD - galaxies: formation - galaxies: kinematics 
and dynamics - galaxies: structure - gravitational lensing 
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1. INTRODUCTION 

Although early-type galaxies appear to be relatively 
simple in terms of their internal structure and stel- 
lar populations, their formation and evolution are still 
poorly understood. In the hierarchical model of struc- 
ture formation, galaxy mergers play a central role in the 
formation and evolution of early-type galaxies. Many 
pieces of evidence are consistent with the merging hy- 
pothesis (see e.g., Ivan Dokkuml 120051 and references 
therein). However, other observational facts appear 
hard to reconcile with substantial amounts of merg- 
ing. Among the facts that have been traditionally 
used to argue against the central role of mergers are 
the high observed central den sity of early-type galaxies 
(|Ostrikerlll980l : ICarlberdll986l ): the small scatter in em- 
pirical scaling laws such as t he color magnitude relation 
([Bower. Lucev fc Ellisl I1992D . the Fundamental Plane 
(here after FP; lDressler et al.lll987t iDiorgovskv fc David 
1987), and the relation between black hole mass Mbh 
and hos t galaxy mass or central stellar velocity op dis- 
persion dFerrarese fc MerrittJl2000tlGebhardt et al.|[2000l : 
IMarconi fc HuntJ|2003t IHaring fc Rixll2004f) : and the cor- 
relation bet ween abundance ratios and host galaxy prop- 
erties (e.g.. lPipino fc Matteucci II2008D . 

Theoretical and ob s ervati o nal studies (e.g., 
IBender. Burstein fc Faberl Il992t ICiotti fc van Albadal 
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2001 1 : iRobertson et al.|[200l lCiotti. Lanzoni fc Volonteril 



2007 : IScarlata et al. r i2007t HFabeT et al . I 120071) have 
argued that a combination of dissipative and dissi- 
pationless processes is needed to match the observed 
properties of early-type galaxies. Dissipation is needed 
for example to create the high density regions and close 
to isothermal density profiles, while dissipationless (or 
"dry") mergers are needed to grow the most massive 
galaxies while preserving the uniformly old stellar 
populations a nd the structural p roperties of present-day 
galaxies (e.g. INaab et all [2007). In particular, both 
theories and observations have suggested that a large 
fraction of the immediate progenitors of today's massive 
red ga laxies could be sp h eroid al galaxies themselves 
( e.g.. Ivan Dokkum et al.l Il999t iK h ochfar fc Burker 



1 20031: INaab Khochfar. fc Burkertll2006t IBell et al.11200 
iFaber et al.l 120071 ). Thus, even tho ugh mergers may not 
be th e only process at work (e.g., IBundv. Treu fc Ellisl 
l2007f) . dry mergers might play an important role, 
especially at relatively recent times, i.e. below redshift 
~ 1.5, after the bulk of star formation is completed for 
massive early-type galaxies in the local Universe (e.g., 
iTrager et al.fl2000l : [Treu et al.|[2005t IRenzini|[200^ . 

In this paper, we take advantage of a new set of con- 
straints provided by strong gravitational lensing stud- 
ies to test the viability of the dry-merger scenario via 
a set of numerical simulations. In practice, we first 
construct realistic models of early-type galaxies so that 
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they obey known scaling relations between their struc- 
tural and kinematic properties, and then study how these 
properties change as a result of dissipationless mergers. 

Our approa ch is similar to tha t fo l- 
lowed bv iNipoti. Londrillo k Ciottil (|2003afl . 
Gonzalez-Garcia fc van Albadal (l2003f) and 

Bovla n-Kolchi n. Ma fc Q uatacrt (2006, see also 
Kazantzidis et al. 2005, Robertson et al. 2006 and 
Johansson, Naab, & Burkert, 2009a for studies that 
include dissipation), who tested how the FP, its projec- 
tions and the A/bh-co relation are affected by mergers. 
The bottom line of these studies is that dry merging 
betw een progen it or sph eroidal galaxies that lie on the 
FP, iKormendvl (|1977t ) and iFaber fc Jackson! (|1976l ) 
relations produce spheroidal galaxies that are close to 
the edge-on FP, but deviate from the Kormendy and 
Faber-Jackson relations, because they have too large 
effective radius and too low central velocity dispersion. 

Taking a step further, we focus here on lensing-derived 
properties, testing the results of our simulations against: 
i) isothermality of the total mass density profile; ii) the 
correlations between lensing mass and size or stellar ve- 
locity dispersion; hi) the Mass Plane, i.e. the correlation 
between surface mass density, effective radius and veloc- 
ity dispersio n discovered by strong len sing studies. As 
discussed bv lBolton et all (|2007ll2008b[ ). the lensing con- 
straints add a qualitatively new piece of information — 
an accurate aperture mass for each galaxy — that helps 
circumvent some of the problems inherent to previous 
studies. The comparison of N-body simulation with tra- 
ditional scaling relations is not straightforward, because 
it relies on assumptions on the mass-to- light ratio to con- 
vert luminosity, which enters the observed correlations, 
into mass. Moreover, among the quantities involved in 
the traditional empirical correlations (e.g., galaxy lumi- 
nosity, size, and velocity dispersion), only velocity disper- 
sion is affected by the galaxy dark matter (hereafter DM) 
content and distribution. Information coming from grav- 
itational lensing, and thus constraining the total galaxy 
mass distribution, effectively solves this problem. 

The paper is organized as follows. First, in §[2] we re- 
view the lensing constraints on the structure of early-type 
galaxies. Then, in §[31 we describe the seed-galaxy mod- 
els, initial orbital configurations, and numerical methods 
that we use to simulate the dry-merging process. In § [4] 
and §[5] we present the results of comparisons between our 
merger simulations and observations in terms of the effect 
of dry mergers on density profiles, central DM fractions 
and scaling relations of galaxies throughout the merging 
hierarchies. In §[B]we discuss our results in the context 
of the redshift evolution of early-type galaxies. Finally, 
§[7] presents our conclusions. 

2. LENSING CONSTRAINTS ON THE STRUCTURE OF 
EARLY-TYPE GALAXIES 

The properties of massive early-type galaxies are 
strongly constrained by the combination of kinematic 
and strong gravitational lensing observations. In partic- 
ular, there is growing evidence that these systems have 
total mass density profile close to isothermal p ex r~ 2 
over a large radial interval (Gavazzi et al. 20071), with an 
intrinsic scatter in logarithmic slope of less th an ~ 10% 
(|Koopmans et aLll2006HKoopmans et al.ir2009f h 

In addition. iBolton et al.l (|2007ll2008bf ) . using a sample 
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Fig. 1. — Mass-space analogs of the Kormendy and Faber-Jackson 
relations as define d by data on the SLA CS early-type gravitational 
lens sample from Bolton et al. (2008a). Squares show the varia- 
tion of half-light radius R c with lensing mass M?„ 1 stars show the 
variation of velocity dispersion cr e 2 with M^ 2 - 
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Fig. 2. — Stellar, dark-matter and total density distributions of 
the seed galaxy (model A; thin curves) and of the end-product 
(thick curves) of the major-merger simulation 2A5po. rm is the 
angle-averaged half-mass radius of the stellar distribution. 

of strong gravitational lenses, have shown that early-type 
galaxies lie on a Mass Plane (MP) 

log R e = a m log er c2 + b m log E o2 + const, (1) 

where a C 2 is the projected velocity dispersion within an 
aperture radius R c /2 and S G 2 is the surface mass density 
within Re/2, with a m = 1.82 ± 0.19, b m = -1.20 ± 0.12 
and RMS orthogonal scatter of 1.24 when normalized by 
the observational errors. The MP is the gravitational 
lensing analogue of the traditional FP. 

The fact that (a m , 6 m ) are consistent with (2,-1) and 
that the scatter is small can be expressed in terms of 
structural and dynamical homology of the lenses, by 
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Fig. 3. — Total density profiles of the seed galaxy model (dotted 
curve) and end-product (solid curve) of a parabolic head-on merg- 
ing between two identical systems with truncated singular isother- 
mal sphere total density distribution, tot is the angle-averaged 
half-mass radius of the total mass distribution. 



defining the dimcnsionlcss structure parameter 



Cc2 



2GM C P 2 



(2) 



where M^ 2 is the total projected mass within R e /2. For 
their sample of lens early-ty pe galaxies from the Sloan 
Lens ACS (SLACS) Survey, IBolton et all (I2008bf) find 
on average 

(log c c2 > = 0.53 ± 0.057. (3) 

We note that the observed scatter on (logc C 2) is 0.08, but 
here we consider the estimated intri nsic scatter . 057 (see 
IBolton et af]|2008bf) . Importantly. IBolton et ail (|2008bf ) 
find that c e 2 has no systematic variation with galaxy 
mass: this is the sense in which the lens galaxies form a 
homologous family. 

In the present paper we also c onsider the analo gues 
of the traditional luminositv-siz e dKormen o^ll977j) and 
luminosity-velocity-dispersion (jFaber fc Jackson! 119761 ) 
relations, by deriving the projected-mass-size and 
the projected-mass-velocity-dispersion relation for the 
SLACS sample of lenses. 

To derive total-mass versions of the Kormendy and 
Faber-Jackson relations, we use the velocity-dispersion, 
half-light radius, and lensing-mass measure ments for the 
SLACS strong-lens sample as published in IBolton et alj 
(20083). We restrict our attention to the 57 lens sys- 
tems wherein the lens galaxy has early-type morphol- 
ogy and single multiplicity. For the mass-velocity dis- 
persion relation, we further restrict ourselves to the 
53 systems for which the Sloan Digital Sky Survey 
(SDSS) spectroscopy is of sufficient sign al-to-noise to 
permit a velocity-dispersion measurement. IBolton et alJ 
(2008a) provide lensing measurements using both sin- 
gular isothermal ellipsoid (SIE) and light-traces-mass 
(LTM) models for the mass density of the foreground 
lensing galaxies, which can be used to gauge the depen- 
dence of results upon the radial mass-density profile used 
to extract the strong-lensing aperture masses. 



Taking the SIE aperture masses as our reference values, 
and fitting for the relation that minimizes the orthogonal 
scatter in logarithmic space, we find a mass-size relation- 
ship given by 



= (iuw ±...n:>) 



x log 



Me P 2 



10 n M, 



o 



0.68 ±0.01. 



(4) 



The observed scatter in logii e at fixed mass is 0.051, 
which is reduced to 0.048 by subt racting in quadrature 
the error in log R as estimated by IBolton et al.l (|2008af ) 
to obtain an estimate of the intrinsic scatter in this quan- 
tity. We note that there will be a degree of covariance 
between the axes of this relation, since the measured R e 
is used to evaluate the aperture-mass of the SIE model 
at R e /2 to obtain M^ 2 . Fitting similarly for the mass- 
velocity-dispersion relation, we find 



log 



0"c2 



100 km s" 



x log 



M c p 2 

10 n M m 



= (0.18 ±0.02) 
-0.36 ±0.01. 



(5) 



This has an observed scatter in logcr o2 at fixed mass of 
0.045, implying an intrinsic scatter of 0.034. Errors in all 
coefficients are computed from the RMS scatter in fits 
to 1000 bootstrap-resampling draws of 57 or 53 galaxies 
apiece. Using LTM lensing aperture masses instead of 
those from SIE models, the coefficients of both relations 
change only insignificantly. These relations are shown in 

Fig.ffl J , 

As discussed i n several papers (Bolt on et al.l 2006: 
iTreu et all 120061 : IBolton et aiTl2008al ITreu et all I2009D 
the SLACS lenses are found to be indistinguishable from 
control samples of SDSS galaxies with the same stel- 
lar velocity dispersion and size, in terms of luminos- 
ity/surface brightness, location on the FP, and distribu- 
tion of environments. We will thus assume that, within 
the range of observational errors, results found for the 
lens sample, including the MP, are generic properties of 
the overall class of massive early-type galaxies. 

3. DISSIPATIONLESS MERGING SIMULATIONS 

We ran N-body simulations of dissipationless galaxy 
merging, exploring both major and minor mergers. In 
the case of major mergers we simulate two colliding sys- 
tems, which are equal-mass, identical galaxy models with 
stellar and DM components, and no gas. In major- 
merger hierarchies, galaxies grow in mass by subsequent 
binary equal-mass mergers, up to an increase in stellar 
mass of a factor of ~ 8. In the first step of a hierarchy the 
seed galaxies are spherically symmetric and non-rotating. 
In the second step of a hierarchy the end-product of the 
first step is replicated, and two such galaxies encounter 
one another. At this stage, the systems are no longer 
spherically symmetric, thus each of them is initially ran- 
domly rotated around its center of mass. The setup of 
the initial conditions of the third step are obtained in the 
same way, starting from the remnant of the second step. 

In the case of minor mergers, we simulate the accretion 
of several smaller satellites onto a more massive system 
(the mass ratio between each satellite and the massive 
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galaxy is either 1/10 or 1/5). For a galaxy to grow sig- 
nificantly on reasonably short timescales through minor 
mergers, more than one satellite is expected to interact 
with the galaxy at the same time. For this reason, we 
considered hierarchies of minor mergers in which each 
step is represented by a multiple merger, so the initial 
conditions represent a massive galaxy surrounded by five 
or ten infalling satellites. Each galaxy collision is fol- 
lowed up to the virialization of the resulting stellar sys- 
tem. A fraction of the particles is unbound at the end of 
the simulation: we define the merging end-product as the 
system composed by the bound stellar and DM particles. 
We note that in the case of minor-merger hierarchies not 
all the satellites present in the initial conditions must 
necessarily merge with the central galaxy. In some cases 
one or two of the satellites are found on an unbound or- 
bit when all the other systems have already merged and 
the merger remnant relaxed. When this happens, the 
satellites are considered escaped and the merging is thus 
characterized by a significant stellar and DM mass loss. 

Altogether we simulated 12 merging hierarchies: 6 
major-merger hierarchies of 3 steps each, 2 major-merger 
hierarchies of 2 steps each, and 4 minor-merger hierar- 
chies of 2 steps each. 

3.1. Seed galaxy models 

The seed galaxy models are spherically symmetric 
equilibrium two-component systems, with stellar and 
DM components. The stellar density dist ribution of the 
seed galaxies is repre sented by a 7 model ()Dehnenlll993l : 
iTremaine et al.lll994[ ): 



P*(r) 



3 — 7 M*f* 
4-7T rT(f* + r) 4 -T 



(0 < 7 < 3), (6) 



where M* is the total stellar mass. We assu me 7 = 1.5, 
thus o ur models are centrally steeper tha n IHernquistl 
(1990) models (7 = 1) and shallower than Uaffel (|1983f ) 
models (7 = 2). We use f* and M* to define 
time and velocity units t* — (f^/GM*) 1 / 2 and = 
f ji* = (GM./f,) 1 / 2 . T h e DM halo is described by 
a iNavarro. Frenk fc White] (|1996l NFW) model, so the 
DM density distribution is 



Pdm{r) 



M, 



dm,0 



r(r + r B )'- 



■ exp 



(7) 



where r s is the scale radius, M^mfi is a reference mass 
and we adopt an exponential cut-off to truncate the dis- 
tribution smoothly at the virial radius r v - n , so the to- 
tal DM mass Mdm = 47r / n °° pd m ( r)r 2 dr is finite. We 
assume Osipkov-Merritt (Osipkovl 119791 : iMerrittl 119851) 
anisotropy in the velocity distribution of the stellar com- 
ponent, whose distribution function is then given by 



f(Q) 



dg* 



V8tt 2 dQ J c^tot VQ - *t 



where 



(8) 
(9) 



The variable Q is defined as Q = £ — L 2 /2r 2 , where the 
relative (positive) energy is given by £ = ty t ot — v 2 /2, v is 



the modulus of the velocity vector, the relative (positive) 
total potential is \& t ot = + ^dm and ^dm are, 
respectively, the relative potentials of the stellar and DM 
components), L is the angular momentum modulus per 
unit mass, and f(Q) = for Q < 0. The quantity r a is 
the so-called anisotropy radius: for r 3> r a the velocity 
dispersion tensor is mainly radially anisotropic, while for 
r <C r a the tensor is nearly isotropic. In the limit r a — > 
00, Q = £ and the velocity dispersion tensor becomes 
globally isotropic. 

The orbital distribution of the DM halo is assumed 
isotropic, so the distribution function of the DM com- 
ponent is given by equation (JSJ) where Q — £ and 
pdm( r ) substitutes f?*(r). The total (stars plus DM) den- 
sity profile is p to t(r) = Pdm(f) + P*i r ) and M tot (r) = 
47r J Q ptot(r')r' 2 dr' is the total mass within r. This fam- 
ily of 7 = 1.5 plus NFW models have four free parame- 
ters: concentration C200 = r V i T /r s , the stellar mass frac- 
tion /» = M„/{M* + M dm ), the ratio £ = r s /R c and 
the anisotropy radius r a . We consider different choices 
of these parameters for the seed galaxies of our hierar- 
chies. In what we refer to as model A (see Table [TJ, we 
assume G200 = 7, /* = 0.09, £ = 7.8, so that the result- 
ing total mass profile is close to isothermal (ptot <x r~ 2 ) 
in the radial range 0.05 <r/rM^ 10 (see thin solid curve 
in Fig. [2|) , where tm is the radius of the sphere contain- 
ing half of the stellar mass, and r a /f* = 1.4, so it lies 
exactly on the MP (logc e 2 — 0.53). We also consider 
models with C200 = 7, /* = 0.09, £ = 7.8, but different 
values of r a /f»: an isotropic model (r a /r* = 00; model 
B) and a more radially anisotropic model (r a /f* = 0.7; 
model C). In other hierarchies the seed galaxy is more 
DM dominated (model D): a 7 = 1.5 plus NFW model 
with C200 = 8, /* = 0.02, £ = 11.6 and r a /f* = 1 so 
that also this model has close-to-isothermal total density 
profile and has logc C 2 ~ 0.53. 

Our seed galaxy models are primarily constructed to 
reproduce the total density distribution and structural 
parameter c C 2 of massive (M* <~ 10 n MQ) local early- 
type galaxies, as constrained by strong gravitational 
lensing and kinematics; we did not attempt a full ex- 
ploration of the space of parameters C200, /* and £, 
but we adopted values that are expected to bracket the 
ran ges expected on the basis of t heoretical investigations 
fsee lNipoti. Treu &; Boltonl2008l and references therein) . 
Additional observational constraints on the DM con- 
tent of present-day early- type galaxies come from weak 
gravitational lensin g surveys (jMandelbaum et ahl 120061 : 
IGavazzi et al.ll2007[). According t o the weak-lensing re- 
sults of iMandelbaum et all (|2006h the stellar mass frac- 
tion /* of M* ~ 10 11 M Q early- type galaxies are of the 
order of 0.03 — 0.06 (if uncertainties on the Initial Mass 
Function are taken into account), so the values of /* of 
our seed galaxy models bracket the observationally esti- 
mated range of /* for systems of this mass. However, /» 
estimated from weak-lensin g tend to decrease wi th stellar 
mass above M* ~ lO n M (jMandelbaum et al.ll20M ). so 
the more DM dominated model D (/» = 0.02) should 
be considered more realistic than model A to represent 
more massive systems. 

To address the question of whether r~ 2 density pro- 
files are preserved by merging we also ran an equal-mass 
merger simulation in which each seed galaxy is a one- 
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component system with total density profile 



Ptot(r) = ptot.o (— ) 



exp 



(10) 



where /9tot,o is a reference density. The density distribu- 
tion in the equation above is that of a singular isothermal 
sphere (SIS), smoothly truncated at a radius rt so that 
the total mass M to t = 47r / °° ptot{r)r 2 dr is finite. This 
truncated SIS model is assumed isotropic. 

3.2. Major-merging simulations 

In the case of major-merging hierarchies, in each step 
we follow the evolution of a binary galaxy encounter, 
characterized by the properties of the colliding systems 
and of their mutual orbit, which we parametrize with the 
standard two-body approximation. For an encounter be- 
tween two galaxies of total masses Mi and M2 we define 
the orbital energy per unit mass 



E = \vl 



G(M 1 + M 2 ) 



(11) 



and the orbital angular momentum per unit mass L = 
do x vo, where do and Vo are the relative separation and 
velocity vectors of the centers of mass of the two collid- 
ing systems at the initial time of the simulation. We ran 
simulations of parabolic (E — 0), elliptic (E < 0) and 
hyperbolic (E > 0) orbits, and we explored both head- 
on encounters (L = 0) and off-center encounters (L > 0). 
In the case of encounters between two equal-mass, iden- 
tical galaxies it is useful to define the dimensionless or- 
bital energy E = 2E/cr 2 and angular momentum modu- 
lus L = L/a v rM.tot, where a v is the virial velocity dis- 
persion and rM.tot is the angle-averaged half-mass ra- 
dius o f the total density distribu tion of the isolated sys- 
tems (|Binnev fc Tremainelll987l ). We consider six differ- 
ent major merging hierarchies, with seed galaxy model 

A, in which the orbital parameters E and L are the same 
in each step of the hierarchy. Only in the last step (run 
8AII10) of the off-center hyperbolic major-merger hier- 
archy E and L are smaller than in the previous steps, 
in order to avoid a too long merging timescale. In addi- 
tion we explore two parabolic hierarchies having as seed 
galaxy model D (one with L = 0, the other with L > 0). 
The parameters of the simulations are summarized in 
Tables [1] and H 

3.3. Minor-merging simulations 

We present four minor-merging hierarchies, each con- 
sisting of two steps of (multiple) minor-merging simula- 
tions. Let us focus first on the hierarchy represented by 
the two simulations 2Alml0x and 4Alml0x: in the first 
step of this hierarchy the initial conditions consist of a 
central galaxy model A, of stellar mass M* and scale 
radius f„, surrounded by ten identical satellites, also 
represented by model A, but with stellar mass M*/10 
and scale radius 0.23f*. The satellites are randomly dis- 
tributed around the central galaxy with relative distance 
between each satellite and the central galaxy in the range 
100 — 140 f * and with velocity directed towards the cen- 
ter of mass of the central galaxy. The relative speed be- 
tween each satellite and the central galaxy is such that 



the orbit would be parabolic in the absence of the other 
satellites. The end-product of run 2Alml0x is used as 
central galaxy of the initial conditions in the second step 
of the hierarchy, in which it is surrounded by ten identi- 
cal satellites of mass M*/5 and scale radius 0.36f* with 
distribution analogous to run 2Alml0x, but now with 
relative distance with respect to the central galaxy in 
the range 200 — 280 r » . Though each satellite would be 
on a radial parabolic orbit in the absence of the other 
satellites, the latter act as perturbers, so the orbits are 
actually neither strictly radial nor strictly parabolic. 

The set-up of the other minor-merging hierarchies 
are similar to the one described above. The hierarchy 
(2Alml0y, 4Alml0y) is just another realization of the 
same set-up: the initial conditions differ just in the rel- 
ative distribution in phase space of the centers of mass 
of the satellites. In the hierarchy (2Alm5, 4Alm5) we 
adopt five satellites per step instead of ten, with stel- 
lar mass 0.2M* (0.4M») and scale radius 0.36f* (0.56f*) 
in the first (second) step. In the hierarchy (2DlmlO, 
4Dlml0) the initial central galaxy and the ten satellites 
are represented by model D, with same scale radii and 
masses as (2Alml0x, 4Alml0x), but relative distance 
between each satellite and the central galaxy in the range 
170 — 230?* (first step) and 340 — 460 f* (second step), 
on account of the more extended DM distributions. 

3.4. Numerical methods 

For the simulations we used the parallel N- 
body code FVFPS (Fortran Version of a Fast 
Poisson Solver; lLondrillo. Nipoti. Ciottil 120031: 
INipoti. Londrillo fc Ciottil l2003al) . based on iDehnenl 
(2002) scheme. We adopt the following values for 
the code parameters: minimum value of the open- 
ing parameter m i n = 0.5 and softening parameter 
e = 0.03 — 0.1 in units of f* (depending on the number 
of particles and on the size of the merging systems). 
The time-step At, which is the same for all particles, 
is allowed to vary adaptively in time as a function of 
the maximum particle density p max : in particular, we 
adopted At = 0.3/(47rGp max ) 1 / 2 . 

Numerically, the initial distribution of the seed galaxy 
particles in phase space is rea lized as described in 
INipoti. Londrillo fc Ciottil (|2002fl . While the gravita- 
tional potential 'i', is analytic, $d m is computed by nu- 
merical integration, given the spherically symmetric den- 
sity distribution ([7]). We verified that the seed galaxy 
models are in equilibrium by evolving them in isolation 
for several dynamical times. As a rule, in the first step 
of each merging hierarchy we use ~ 4 x 10 4 stel- 
lar particles and A dm ~ 4 x 10 5 (when /* = 0.09) 
or iV dm ~ 2 x 10 6 (when /» = 0.02; halo and stel- 
lar particles have the same mass). Clearly, the num- 
ber of particles increases with the hierarchy steps, up to 
N tot = N dm + ~ 1.7 - 3.9 x 10 6 in the last steps 
considered. 

The intrinsic and projected properties of 
the merger remnants are d eterm ined following 
INipoti. Londrillo fc Ciottil (|2002l l2006f ). In particu- 
lar, we measure the axis ratios c/a and b/a of the 
inertia ellipsoid (where a, b and c are the major, 
intermediate and minor axis), and the angle-averaged 
density distribution of the stellar and DM components. 
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We also measure the angle-averaged half-mass radius 
of the stellar density distribution (Ym) and of the total 
density distribution (rM.tot)- For each end-product, in 
order to estimate the importance of projection effects, 
we consider 50 random projections and we use them to 
compute average projected quantities and associated 1-cr 
scatter. For each projection we measure the ellipticity 
e = 1 — b c /a c , the circularized projected density profile 
and the circularized effective radius R c = yja c b c (where 
a c and b c are the major and minor semi-axis of the 
effective isodensity ellipse). In addition, we measure 
the central projected velocity dispersion cr c2 , obtained 
by averaging the projected velocity dispersion over the 
circularized surface density profile within an aperture of 
Re/2. The projected mass M£ 2 is computed by counting 
stellar and DM particles within the same projected 
elliptical aperture with circularized radius R e /2. 

As a convergence test, we reran three of the simulations 
at higher resolution (with five times more particles), find- 
ing good agreement with the lower-resolution cases. In 
particular, the runs 2A5ph, 2A5po and 2A5ml0x differ 
from the runs 2Alph, 2Alpo and 2Alml0x, respectively, 
just for the number of particles of the seed galaxies (see 
Tables Q] and \2§ : the quantities R e , a e2 and M^ 2 , which 
are crucial for the present study, are indistinguishable in 
lower and higher resolution simulations. 

4. EFFECT OF DRY MERGERS ON THE MASS DENSITY 
DISTRIBUTION 

4.1. Total density profiles 

From lensing studies we know that the total density 
profiles of early-type galaxies are close to isothermal over 
a large radial range. Here we address the question of 
whether such a feature is preserved by dry major and 
minor merging. As a first experiment we consider the 
merging between two identical truncated SISs with total 
density distribution (TTU)) . In this simple case each seed 
galaxy model is one-component, so we do not distinguish 
between DM and baryons. Figure [3] plots the total den- 
sity profiles of the seed models and of the merger remnant 
of this run. By construction, the seed model has den- 
sity profile oc r~ 2 within ~ J~M,tot and steeper at larger 
radii because of the exponential truncation. The den- 
sity profile of the merging end-product is very similar to 
that of the seed galaxy: fitting it with a power law r -7 
in the radial range 0.03 < r/rM,tot < 1 we get a best- 
fitting index 7' = 2.00 ±0.01. The finding that the inner 
density cusp is preserved by dissipationl ess merging is 
consistent with theoretical expectations (Dchncn 2005) 
and previous numerical experiments (IFulton fc Barnea 
I200U iKazantzidis. Zentner fc Kravtsovll2006ft . Our sim- 
ulation shows that the ~ r~ 2 density profile is preserved 
by dry major merging not only in the central regions, but 
throughout the system. 

Let us now consider how the total density profile is af- 
fected by merging more realistic (two-component) stel- 
lar systems. Figure 0] plots, as a function of radius, 
the (angle-averaged) total density distributions of our 
six reference major merging hierarchies, having as seed 
galaxy model A and differing for the orbital parameters 
of the galaxy encounters. From the diagrams it is ap- 
parent that the merging process affects significantly the 
normalization of the density profile (producing systems 



of lower density at the half-light radius; see also § [5]), 
but it roughly preserves its shape. To quantify the ef- 
fect of dry merging on the shape of the total mass den- 
sity distribution we fit the total intrinsic density pro- 
files of our models with a power law r -7 over the ra- 
dial range 0.1<r/i? c <l, where R c is the average of 
the circularized half-light radius over 50 random pro- 
jections. By construction our seed galaxy models have 
close-to-isothermal total density profile: in the radial 
range 0.1 < r/R c < 1 model A and model D are best fitted 
with 7' = 2.02 ± 0.02 and 7' = 2.00 ± 0.01, respectively. 
Figure [5] shows the distribution of the best- fitting 7' for 
the end-products of all steps of our merging hierarchies 
(30 steps altogether). It is apparent that dry merging 
can make slightly steeper or shallower the total density 
profile, depending on the characteristics of the galaxy 
encounters. The distribution has mean (7') ~ 2.01 
and standard deviation oy ~ 0.18 (~ 9%). This is 
comparable with the observationally determined values 
7') ~ 2.08 ± 0.1 and standard deviation <r 7 < < 0.20 
Koopmans et al.lf2009f h suggesting that if a population 
of early-type galaxies has initially close-to-isothermal to- 
tal density profile, dry merging tends to maintain this 
property, but adds non-negligible scatter. 

The best-fitting values of 7' for the single simulations 
are reported in Table [2] together with the other simula- 
tion parameters: off-axis major mergers tend to produce 
steeper density profiles than head-on encounters, while 
there is no evidence of a dependence of 7' on the orbital 
energy. Minor and major merging involving more DM 
dominated galaxies (model D) tend to decrease 7'. 

We focused on the slope of the total density profile 
over the radial range 0.1<r/i? c <l, which is approxi- 
mately the range probed in combined strong-lensing and 
stellar kinematics observations. However, we note that 
the total density profiles of our seed models and merger 
remnants are well represented by a power-law over a 
much larger radial range, though the power-law index is 
rather sensitive to the radial interval. For instance, seed 
models A and D are fitted over the larger radial range 1 
0.1 <r/R c < 10 with 7' = 2.17±0.02and7' = 1.94±0.02, 
respectively. Also fitting over this larger radial range, the 
end-products are found to have slightly larger or smaller 
best-fitting 7' than the seed galaxy of their hierarchy. 

4.2. Central dark-matter fractions 

We have seen that dry merging tends to preserve the 
slope of the total density profile. We also found that the 
shape of the density distributions of stars and DM is not 
dramatically affected by dry mergers (see Fig. , though 
there is redistribution between DM and stars that leads 
to a systematic increase with merging of the projected 
DM-to-total mass ratio /? within R c . This is apparent 
from the upper panel Fig. [&1 plotting f% m (R < R e ) as a 
function of M* for our major and minor merging hierar- 
chies. The seed galaxy models have /? = 0.25 (model 
A) and /? = 0.45 (model D), while the merger rem- 
nants have DM fractions in the range 0.3</j <0.7. 
In all cases increases with the stellar mass along 
the merging hierarchy. Taken by itself this effect might 

1 Larger radial ranges are interesting for comparison with we ak- 
lensing observations of early-type galaxies (Gavazzi et al. 2007). 
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Fig. 4. — Total density distributions of the seed galaxy (solid curves) and of the end-products of six major-merging hierarchies having 
the same seed galaxy (model A), but different orbital parameters. Re is the average of the effective radius over 50 random lines of sight. 
The end-products of the first, second and third steps of the hierarchies are represented by dotted, short-dashed and long-dashed curves, 
respectively. The dash-dotted line represents the isothermal distribution (ptot on r~ 2 ) in the radial range 0.1 < r/R c < 1, over which we 
fit the profiles. 



be interpreted as quite naturally explaining the "tilt" 
seen in the luminosity-based FP relation in terms of a 
systematic increase in the central DM fraction with in- 
creasing mass. If more luminous and massive galaxies 
represent a later stage in the merging hierarchy, this 
increased central DM fraction will lead to a larger dy- 
namical mass per unit luminosity as compared to less 
massive (and less merged) galaxies. However, the in- 
crease of f^ m {R < R c ) with mass during the merging 
hierarchy is mainly driven by a strong increase in the 
effective radius produced by dry merging, which might 
be considered an argument against the relevance of dry 
mergers in the assembly of elliptical galaxies (see § [5]). 
In fact, the trend is quite different if one considers 
as measured within a small aperture of fixed radius, in- 
dependent of the galaxy's half-light radius. For instance, 
in the lower panel of Fig. [6] we show the results ob- 



tained by measuring the projected DM fraction within 
and aperture of radius R = r*. We note that f» ranges 
from ~ 0.8R e for the seed galaxies, to f* ~ 0.07i? o for 
the most extended remnant. We find that the remnants 
have fd m {R < r*) close the value measured in the seed 
galaxy of their hierarchy, indicating that the central pro- 
jected DM fraction, measured within a region of fixed 
size in physical units, is not substantially affected by dry 
merging. Interestingly, minor mergers (stars, diamonds 
and crosses in the diagrams) tend to produce systems 
with higher values of f^ m than major mergers (circles 
and squares in th e diagrams), indicat i ng that dynamical- 
friction heating (lEl-Zant et alj 12004 INipoti et alj|2004t 
iRomano-Diaz et al.l 20081 Uohansson. Naab. fc Ostriken 
2009b) is effectivel y balanced by accre tion of DM in the 
central regions (see INipoti et aLlfeOO-H . for a discussion). 
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Fig. 5. — Distribution of the values of the best-fit power-law index 
of the total density distribution ptot oc r ' for the end-products of 
all steps of minor and major merging hierarchies. The distribution 
has mean {7') ~ 2.01 and standard deviation ay ~ 0.18. 
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Fig. 6. — Upper panel: projected dark-matter-to-total-matter 
mass ratio within ftasa function of the stellar mass for the major 
(circles and squares) and minor (stars, diamonds and crosses) merg- 
ing hierarchies. The points representing the seed (M*/M„ = 1) 
galaxy models A and D are labelled. For the merger remnants 
depends on the line of sight, but projection effects are small, 
so just the average value is plotted. Lower panel: same as upper 
panel, but plotting projected dark-matter-to-total-matter mass ra- 
tio within an aperture of radius Ft = f * (same radius, in physical 
units, for all models). 

5. EFFECT OF DRY MERGERS ON THE LENSING 
SCALING RELATIONS 

We address here the question of whether the lens- 
ing scaling relations are preserved by major and mi- 
nor dry merging. For this reason we consider seed 
galaxy models that lie on the observed correlations and 
we investigate whether the merger remnants are consis- 
tent witdT_the_sjmT£j^ intrinsic scat- 
ter. iNipoti. Londrillo & Ciottil (|2003af) pointed out that, 
when comparing merging simulations with the standard 



(luminosity based) scaling relations of early-type galax- 
ies, it is crucial to consider not only the edge-on FP, but 
also the Faber- Jackson and the Kormendy relation, be- 
cause a system may well lie on the FP while deviating 
significantly from the other relations. Similarly, the fact 
that a galaxy model lies on the MP does not mean that it 
satisfies also the the observed lensing mass-size and mass- 
velocity dispersion relations, because — for fixed M p 2 — 
anomalous values of both cr e 2 and R e can conspire to 
give a value of c e 2 perfectly consistent with observations. 
For this reason we will study how our merger hierarchies 
behave with respect to all three lensing scaling laws. 

5.1. Placing the seed galaxies on the observed scaling 
relations 



As discussed extensively in INipoti. Treu &: Boltonl 
(2008), for a spherical galaxy model the dimensionless 
quantity c e 2 is fully determined by the stellar and DM 
density profiles, and by the distribution of the orbital 
velocity of stars. Of course, purely dynamical systems 
such as our seed galaxy models can be rescaled arbitrar- 
ily by choosing their mass scale M* and length scale f * in 
physical units. While the value of c C 2 is independent of 
this choice, in order to place a model on the observed 
mass-size and mass-velocity dispersion relations, mass 
and length scales must be chosen consistently. 

Given that the MP relation has no tilt and that the ob- 
served mass-size and mass-velocity dispersion relations 
are consistent with power laws over the range covered 
by the lensing data (equations [Tj , [3] and [5] ) , we would 
not need to specify the mass of the seed galaxy in or- 
der to address the question of the effect of dry merg- 
ing. It is sufficient to check whether the merging end- 
products have values c e 2 consistent with the MP and the 
relative change of M^ 2 , R c and er c2 are consistent with 
the observed mass-size and mass-velocity dispersion re- 
lations. However, just for clarity and without loss of 
generality, we fix M* = lO n M0, so all our seed galaxies 
have total stellar mass M* = 10 11 M Q . The projected 
mass M c p 2 = 3.73 x 1O 1O M for galaxy model A and 
M c p 2 = 4.48 x 10 10 M Q for galaxy model D. We then fix 
f* = 2.08 kpc for model A and f* = 2.34 kpc for model 
D, so that they lie on the M p 2 — R e relation (and, as a 
consequence, on the M p 2 — er C 2 relation, because they lie 
on the MP by construction). The satellites in the minor- 
merger hierarchies have masses and characteristic radii 
proportionally smaller than those of the central galaxy 
(see § 13. 3[) . but such that they also satisfy the scaling 
laws. Clearly, M* and f* are kept fixed throughout each 
hierarchy, so masses, sizes and velocities of all remnants 
are now known in physical units. 

5.2. Comparing the merger remnants with the observed 
scaling relations 

The location of the merger remnants with respect to 
the observed scaling laws is shown in Figs. [3 [8] and [9l 
We note that the total projected mass M p 2 within R e /2 
increases with dry merging more than the stellar mass 
or the total mass integrated over the entire system. This 
can be seen, for instance, in Fig. [7] in three merging 
steps M p 2 increases by more than a factor of ten, while 
the stellar or total mass, integrated over the entire sys- 
tem, can increase at most by a factor of eight (in fact 
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Fig. 7. — From top to bottom: best-fitting logarithmic slope 7' (pt. 



) of the total density profile (first row), structure parameter 



c e 2 (second row), effective radius (third row) and projected velocity dispersion <r e 2 (fourth row), as functions of the projected mass M c p 2 , 
for elliptic (left column), parabolic (central column) and hyperbolic (right column) major-merging hierarchies with seed galaxy model A. 
Points and error bars represent the average value and the i-tr scatter (due to projection effects) for simulated galaxies: filled circles are for 
head-on encounters and empty circles for off-axis encounters. M^„ depends on the line-of-sight, but error bars are smaller than the symbol 
size, so they are not plotted. Solid and dotted lines represent the observed correlations with 1-a scatter. 



less, as a consequence of mass loss). This finding indi- 
cates that dry merging redistributes matter so that the 
fraction of mass found in a cylinder of radius R e /2 in- 
creases significantly. 

Let us focus now on the behavior of the merger hier- 
archies with respect to the MP, which is illustrated in 
the second row of panels of Fig. \7\ [5] and In general, 
the intermediate and final remnants of all our merger 
hierarchies are found to be very close to the MP. Ac- 
tually, all the remnants of model-A merging hierarchies 
have values of c C 2 within the observed scatter of the MP, 
while the end-products of the model-D merger hierarchies 
(squares and crosses in Fig. [S]) have slightly higher values 
of c C 2 , and are only marginally consistent with the MP. 
This result implies that c C 2 is a robust parameter, since 
its value is essentially preserved under a wide range of 



dry-merging histories. In a scenario in which early-type 
galaxies grow significantly by dry merging, the charac- 
teristic value of c e 2 must therefore be established by the 
dissipative ( "wet" ) formation processes of the primordial 
progenitor galaxies. The implications of the observed 
value of c C 2 for galaxy s tructure (as opposed to evolu - 
tion) are the subject of iNipoti, Treu fc Boltonl (|2008l) . 
Remarkably, dry merging appears also to wash out the 
effect of orbital anisotropy on c e 2'. the end-products of 
mergers whose seed galaxies differ only in the velocity 
distribution of their stars (and thus ha ve different values 
of c 2; see lNipoti. Treu &: Boltonll2008f ) are indistinguish- 
able in the space of the parameters M^ 2 , R e and cr e 2 (see 
Fig.®. 

As anticipated above, the fact the MP is "closed" 
with respect to dry merging does not mean that 
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Fig. 8. — Same as Fig. [7] but for parabolic major-merger hierar- 
chies having as seed galaxy model D (solid squares are for head-on 
encounters and empty squares for off-axis encounters), and for the 
following minor merging hierarchies: model-A hierarchies with ten 
satellites per step (stars) and five satellites per step (diamonds), 
and the model-D hierarchy with ten satellites per step (crosses). 

the same must be true also for the mass-velocity- 
dispersion and mass-size relations. In fact, as 
might be expected from previous explorations of 
the effect of dry mergin g on the Faber- Jackson 
and Kormendy relations (Nipoti. Londrillo k. Ciottil 



I2003al iBovlan-Kolchin. Ma k Quataertl 120061) . the ench 
products of all the merger hierarchies are found to devi- 
ate significantly from the M^ 2 -R c and M^ 2 -a c2 relations, 
because they have too large R c and too low cr C 2 for their 
mass (see third and fourth rows of panels in Fig. [7j [8] 
and [5]). It turns out that the too large R e and too low 
<7 C 2 compensate, giving values of c e 2 very close to the ob- 
served one (we recall that R e and a e 2 appear in the defi- 
nition of c e 2 just in the product R c a^ 2 \ see equation [5]). 
Some of the minor-merger hierarchies lead to strong devi- 
ations from the observed scaling relations (see right-most 
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Fig. 9. — Same as Fig. [7] but for the seed models and end- 
products of three head-on parabolic major-merging simulations dif- 
fering only for the anisotropy of their seed galaxy models: 2A5ph 
(circles; r a /r* = 1.4), 2B5ph (triangles; r a /f* = 00) and 2C5ph 
(pentagons, r a /f* = 0.7). 

stars in Fig. [5]): in some steps R c increases with M£ 2 
more than linearly, and a e 2 decreases with Mf 9 , consis- 
tent w ith recent results bv lNaab. Johansson, fc Ostrikerl 
(2009). In other cases the increase in size is significantly 
smaller, suggesting that the structural effect of minor 
merging is quite dependent on the details of the interac- 
tion, such as DM distribution and orbits of the satellites. 
In the merging hierarchies having as seed galaxy model 
D, whose end-products are found slightly off the MP, the 
deviations from the mass-size and mass-velocity disper- 
sion relation are in general relatively small (squares and 
crosses in Fig. [8]), and marginally consistent with the 
observed relations in the case of off-axis major mergers 
(open squares in Fig. [5]). In these cases low values of a e 2 
are not compensated by large enough R c . 
In the top panels of Fig. [7j [8] and [9] the total-density 



Dry mergers of early-type galaxies 



11 




10.5 11 11.5 

log M|,/M e 



Fig. 10. — Top: the green shaded area represents the expected 
distribution, within 1 cr, of the merger remnants of a population 
of galaxies all with the same projected mass M e p 2 = 1O 1O - 6 M , 
and a distribution of c e 2 matching the observed one (see §0. The 
red shaded area represents the location in the M^ 2 — c e 2 plane of 
possible dry-merging progenitors of a z = early-type galaxy with 
projected mass M^ 2 = 10 116 Mq (see §|6}- Central and bottom 
panels: same as top panel, but in the M^ 2 — and M? 2 — °e2 
planes. Symbols are the same as in Figs.[7]and[8] Solid and dotted 
lines represent the observed correlations with l-cr scatter. 
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Fig. 11. — Effective radius (upper panel) and projected velocity 
dispersion (lower panel) of the seed galaxies and merger remnants 
of all steps of our major and minor merging hierarchies. For each 
hierarchy, iff* and f * are set so that the remnant of the last step 
lies on the local iW?„ — R c relation. Solid and dotted black lines 
represent the local observed correlations with 1-<t scatter. Symbols 
are the same as in Fig. [7] and [8] 

slope 7' (see § 01 is plotted as a function of mass for the 
seed galaxies and the remnants. The behavior on the 
M^ 2 -j' plane is not strictly correlated to that in the M^ 2 - 
c C 2, M^ 2 -R e and M^ 2 -a C 2 planes, though we note that 
models that deviate more from the M^ 2 -R e and M^ 2 -a C 2 
correlations tend to have steeper total density profiles 
(i.e. higher values of 7'), while systems with higher val- 



ues of c C 2 have shallower total density profiles. 

5.3. Mass dependence of R c , o~ C 2 and c C 2 for merger 
remnants 

The growth of galaxies by dry merging of smaller 
systems can be quantified by approximating the de- 
pendence of the structural parameters on mass with 
power-law functions: R c cx (M^ 2 ) aR , cr c2 oc {M^ 2 ) a " 
and c C 2 cx (M^)"'. For each step of a hierarchy the 
power-law indices are given by an = A log R e /A log , 
a a ee A log cr o2 / A log M e p 2 and a c = A log c o2 / A log Alf 2 , 
where A indicates variation between the remnant and 
the progenitor [e.g. Alogi? c = (log J R G ) romnant - 
(log Re) progenitor]- We can obtain an estimate of the 
expected global effect of a general dry-merging hierar- 
chy by considering the average values ((or) = 0.85, 
(a a ) = 0.06, (a c ) — 0.04) and the associated standard 
deviations (Soir — 0.17, Sa a = 0.08, Sa c = 0.07) of these 
values for our major and minor merging hierarchies (30 
steps altogether). These values must be compared with 
the slope of the observed scaling laws (equations [3], [4] 
and [S]), which correspond to aji = 0.63, a a = 0.18 
and a c = 0, consistent with our statement that the MP 
is roughly preserved by dry mergers, but the mass-size 
and mass-velocity-dispersion relations are not. 

Remarkably, the distributions of an, a a and a c for our 
N-body experiments are characterized by a substantial 
scatter, which is quantified by the quoted values of the 
standard deviations. Such a broad distribution around 
the average slope is an additional problem for the dry- 
merging scenario, because it implies that dry mergers bf 
would introduce a lot of dispersion in the scaling rela- 
tions. To quantify this, let us consider a population of 
galaxies, all with M e p 2 = Mq and R c = Rq, that evolve 
through dry merging with an unspecified merging his- 
tory made of both minor and major mergers (we are 
not assuming that these seed galaxies merge just among 
themselves): under the hypothesis that an is distributed 
normally, the merger remnants of mass M^ 2 > Mq are 
expected to have logi? e normally distributed with aver- 
age (log -R e ) = logi?o + l°g(M?2/-^o) and standard 
deviation dan\og(M^ 2 / A1q). Thus, in this simple exam- 
ple, the predicted scatter in the M£ 2 ~ Rc relation at a 
given M^ 2 increases for increasing 5a and — for fixed 
San — increases for increasing M e p 2 , suggesting that a 
large enough mass growth by dry merging would pro- 
duce a spread in R c at given M^ 2 inconsistent with the 
tight observed correlation. The importance of this ef- 
fect can be appreciated by looking at the green shaded 
area in the central panel of Fig. [TD1 corresponding to the 
expected distribution, within 1 cr, of the merger rem- 
nants of a population of galaxies all with the same pro- 
jected mass M^ 2 — 10 10 ' 6 M Q (green vertical bar in the 
diagram), but now with a more realistic distribution of 
R c : we assume that logi? e is normally distributed with 
average (log R c ) satisfying the observed M^ 2 -R e correla- 
tion (equation [4]) and standard deviation equal to its 
intrinsic scatter (0.048 dex). Similar considerations ap- 
ply to c C 2 and cr 2 (green shaded areas in the top and 
bottom panels of Fig. [T0|) . It is apparent that the scatter 
in R c and er C 2 for the remnant population rapidly in- 
creases with mass above the observed scatter, while the 
effect is weaker in the M c p 2 -c e 2 plane. The exact shapes 
of the green shaded areas depend on the set of merging 
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hierarchies considered. For instance, in the M^ 2 -R e , we 
find a weaker deviation from the observed correlation if 
we consider only model-D hierarchies, and a stronger de- 
viation if we consider only minor-merging hierarchies (as 
a consequence of the specific behavior of these models; 
see § I5.2p : in any case regions of the diagrams out of the 
observed strip are quickly populated. 

Summarizing, we showed that the lensing mass-size 
and mass-velocity-dispersion relations are effectively de- 
stroyed by dry mergers, which produce systems with too 
large half-light radii and too low velocity dispersions, and 
introduce substantial scatter in R e and cr e 2 at given mass. 
We conclude that present-day early-type galaxies did not 
form by dry merging of high-z early-type galaxies obeying 
the local (z = 0) lensing scaling relations. 

6. DRY MERGERS AND REDSHIFT EVOLUTION OF 
EARLY-TYPE GALAXIES 

In the previous section we explored the effect of 
dry mergers on galaxies lying on the lensing scaling 
relations observed locally. However, higher-z early-type 
galaxies must not necessarily obey the z — lens- 
ing scaling laws. In fact, though the lensing scaling 
relations at higher redshift are not available, mea- 
sures of the effective radii and the stellar masses of 
high-redshift (z > 1 — 2) early-type galaxies suggest 
that these objects may be remarka bly more compact 
than their lo c al counterparts (e . g.. iDaddi et al.1 120051: 
Truiillo et all l2006t IZirm et alj 120071: iCimatti et al 



20081 Ivan Dokkum et all 120081 : Ivan der Wei et al 



20081 ) . As shown by our current results, as well 



as previous work (iNipoti. Londrillo fc Ciottil l2003at 
iBovlan-Kolchin. Ma fc Quataertll2006l ). dry merging has 
the effect of making galaxies less compact, so it has been 
proposed that evolution of early-type galaxies through 
such a mechanism might be consistent with the recent 
observational fi nding that highe r -z ga la xies are more 
compact (e.g.. iKhochfar fc Silk) [2001 iHopkins et all 
12009b Ivan der Wei et al.l l2009t l Even allowing for the 
possibility that the candidate dry-merging progenitors 
of present-day early-type galaxies do not obey the local 
lensing scaling laws, the tightness of these correlations 
can be used to constrain the contribution of dry mergers 
in the assembly history of early-type galaxies. 

For this purpose it is useful to consider the follow- 
ing question: under the hypothesis that local early-type 
galaxies (satisfying the lensing scaling laws) formed by 
dry merging, what are the expected properties of their 
(high-redshift) progenitors? To address this question we 
assume that the remnants of the last steps of all our 
simulated merging hierarchies satisfy the local M^ 2 — R c 
relation. We then fix M* = IO^Mq, and we choose a dif- 
ferent value of f* for each hierarchy (in all cases smaller 
than the values of f* adopted in §[2)), so M e p 2 , R e and a e2 
are given in physical units for all models. The location 
of the models with respect to the observed lensing cor- 
relations is shown in Fig. 111! the most massive systems 
lie (by construction) on the observed M£ 2 — R c relation 
(upper panel) and approximately also on the M^ 2 — o"c2 
relation (lower panel), consistent with the fact that c e 2 
is not significantly affected by dry merging; less mas- 
sive systems lie systematically below the M^ 2 — R e re- 
lation and above the M^ 2 — cr e 2 relation. In this pic- 
ture the lowest-mass systems in the diagrams in Fig. [TTJ 



(those with log M^/M© < 10.7, all having by construc- 
tion M* = 1O 11 M ) would represent the high-z progen- 
itors of local ellipticals. Given the idealized nature of 
our experiments, which are not set in a full cosmological 
context, we cannot define precisely the redshift of these 
progenitors. However, a lower limit to their redshift can 
be obtained by noting that the time to complete our 
merging hierarchies is in the range 6—12 Gyr, with me- 
dian ~ 9 Gyr. A realistic merging hierarchy is expected 
to be a mixture of steps of the explored hierarchies, so 
we take as reference time the median, which is the look- 
back time to z ~ 1.4. Therefore, in the present context, 
our seed galaxy models should be compared with z> 1.4 
early-type galaxies with stellar mass ~ IO^Mq. The di- 
agrams in Fig. [TTJ show that the candidate progenitors 
are predicted to be compact, with effective radii up to 
a factor of ~ 2.5 smaller and velocity dispersions up to 
a factor of ~ 1.6 higher than local galaxies with similar 
mass, and are characterized by a significant spread in 
both R e and a C 2- We recall that the candidate progen- 
itors have values of c e 2 consistent with those observed 
in present-day galaxies, because our seed galaxies lie on 
the MP by construction. Though one cannot assume a 
priori that high-z progenitors lie on the MP, this assump- 
tion is justified as far as we are assuming — as working 
hypothesis — that local early-type galaxies formed by dry 
merging, because we have shown that dry mergers tend 
to preserve the MP. 

The average properties of the candidate progenitors 
are qualitatively c onsistent with those of z<*, 1.4 early- 
type galaxies (see ICimatti et al.l 120081 : iCappellari et al.l 
2009), so one might be tempted to conclude that a dry- 
merging scenario is compatible with observations. This 
is not necessarily the case, because the tightness of the 
local relations must be reproduced as well. In fact, as 
we are going to show, the tightness of the local relations 
sets the strongest constraints on dry merging, requiring a 
remarkable degree of fine tuning of the mix of progenitors 
and types of interaction. 

This can be seen from the example shown in Fig. 1101 
where we focus, without loss of generality, on local early- 
type galaxies with projected mass M e p 2 = 1O 116 M0 (red 
vertical bars in the diagrams), and we constrain the prop- 
erties of their progenitors in a dry-merging scenario. We 
have seen in § 15.31 that progenitors of given M^ 2 , o~ C 2 
and R c would produce through dry merging a popula- 
tion of galaxies with a range of values of cr e 2, Re and 
c e 2 for given mass, depending on the expected distribu- 
tions of ttfl, a a and ot c , which are constrained by our 
N-body simulations. Here we take the opposite approach 
and impose that our reference present-day galaxies with 
M^ 2 = 10 116 Mq have logc C 2, logi? e and loga e 2 with 
average values and intrinsic scatter consistent with ob- 
servations. This procedure delivers constraints on M^ 2 , 
er C 2, -Rc and c G 2 in the form of "allowed" regions in the 
diagrams of Fig. [TU] (red shaded areas). For instance, 
let us focus on the M^ 2 — R c plane (central panel): at a 
given mass M^ 2 < 1O 116 M , the upper and lower limits 
of the red shaded areas bracket the distribution within la 
of log R c of a normally distributed population of galaxies 
of mass M^ 2 that can produce via dry merging our popu- 
lation of galaxies with lensing mass 1O 116 M . The same 
considerations apply to the M<? 2 -c C 2 and M ( f 2 -cr G 2 planes. 
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It is clear that the smaller the shaded area, the higher the 
fine tuning necessary to obtain the local distributions. 

While the constraints in the c e 2 plane are weak (be- 
cause c e 2 is roughly preserved by dry merging), those in 
the <7 e 2 and R c plane are quite strong and indicate that 
an uncomfortable degree of fine tuning is required in a 
dry merging scenario. Note, in particular, that a growth 
in M^ 2 via dry merging by more than a factor of ~ 2 
(and in size by more than a factor of ~ 1.8) is excluded 
(see the red shaded area in the M^ 2 -R c plane). Though 
our results constrain the projected mass M^ 2 measured 
within i? c /2, the upper limits on the mass growth can 
be extended to the total mass, because in our simula- 
tions Mg 2 increases with dry merging faster than the to- 
tal mass (integrated over the entire system; see § 15. 2|) . 
Therefore, our results indicate that present-day, massive 
early-type galaxies did not assembled more than 50% of 
their total (dark plus luminous) mass through dry merg- 
ing. 

One important caveat is that the allowed amount of 
mass growth depends on (or) and San, and thus on the 
considered merging histories. Our set of simulations is 
biased in favour of major mergers and our exploration of 
the orbital parameter space is limited, especially in the 
case of the more DM dominated (model D, /* = 0.02) 
galaxies. For instance, if we consider only our minor- 
merger hierarchies the upper limit to the mass fraction 
goes down to ~ 40%, while this figure is ~ 60% if we 
consider only model-D hierarchies. In any case, a certain 
amount of fine tuning is required. Larger sets of numer- 
ical experiments, based on more realistic merging histo- 
ries, are needed to quantify whether the amount of fine 
tuning required is consistent with hierarchical growth in 
a cosmological context. 

From Fig. llOl it is also apparent that even an increase in 
mass by a factor smaller than 2 requires a lot of fine tun- 
ing, in the sense that the progenitors must have tighter 
mass-size and mass-velocity-dispersion correlations than 
local galaxies. This cannot be excluded, because there 
is no information so far on the lensing scaling relations 
beyond the redshift range covered by the SLACS sample 
(z~ 0.2 -0.3). 

7. CONCLUSIONS 

We explored the effect of minor and major dry mergers 
on the lensing scaling relations of early- type galaxies. 
Our main findings can be summarized as follows. Major 
and minor dry mergers: 

1. Preserve the nearly isothermal structure (ptot oc 
r~ 2 ) of early- type galaxies, adding a scatter of 
~ 9% to the logarithmic slope (for an increase in 
stellar mass of up to a factor of ~ 8), consistent 
with the observed value of ~10%. 

2. Do not change the ratio between total (lensing) 
mass M^ 2 and "virial" mass (R c a^ 2 /2G) more than 
the observed scatter of 0.057 dex. In other words, 
they move galaxies along the MP. 

3. Move galaxies away from the observed correlations 
between total (lensing) mass and size or velocity 
dispersion. Specifically, dry mergers increase the 
radius more rapidly than the scaling relation with 



mass would predict, while they do not increase ve- 
locity dispersion rapidly enough. 

4. Add substantial scatter in the lensing mass-size and 
mass-velocity-dispersion relations . 

The first two findings indicate that two important 
properties of early-type galaxies — the isothermal total 
mass profile and the MP — are quite robust against dry 
merging. If these regular structural and kinematic prop- 
erties are established — presumably at z > 1, in order to 
satisfy con straints on the ol d age of their stellar popula- 
tions (e.g. jTreu et al.ll2005l ) and of the lack of evolution 
of the logarithmic slope of the total mass den sity profile 
(jKoopmans et a l. 2006; Koopma ns et al.ll2009l ) — further 
growth via dry merging would not spoil this structural 
and kinematic homology of early-type galaxies. 

The third and fourth findings pose a significant prob- 
lem for a scenario where early-type galaxies grow signifi- 
cantly in mass via dry mergers. The non-conservation 
of the mass-size and mass-velocity-dispersion relation 
shows that the progenitors of present day early-type 
galaxies cannot be obeying the same scaling relations. 
On average, they will have to be smaller in size 
for a given velocity dispersion or total mass, quali- 
tatively in line with som e recent observational results 
(jvan der Wei et all l2008h . However, and most impor- 
tantly, dry mergers add substantial scatter to the scal- 
ing relations. A few random generations of dry merg- 
ers would be sufficient to produce a scatter much larger 
than observed. In particular, our results exclude that 
present-day massive early-type galaxies assembled more 
than ~ 50% of their total (luminous plus dark) mass and 
increased their size by more than a factor ~ 1.8 via dry 
merging. Quantitatively, the amount of dispersion intro- 
duced in the scaling laws by dry merging depends on the 
detailed properties of the merging history, such as DM 
content of the progenitors and mass-ratio and orbital pa- 
rameters of the galaxy encounters. A high degree of fine 
tuning - where the properties of the progenitor correlate 
with the peculiarities of the mergers process - would be 
needed to reconcile dry mergers with the tight observed 
lensing scaling relations. 

In the present work we considered only gas-free merg- 
ers, so our results do not apply if mergers occur in 
the presence of significant amounts of gas. It is well 
known that dissipative processes might help recon- 
cile merging with the standard scaling relations (e.g . 
Robertson et alj|2006t iCiotti, Lanzoni fc Volonter1 l2007l: 
Johansson. Naab. fe Burkertl |2009eJT " so it is reasonable 
that the same can happen for the lensing scaling rela- 
tions. However, the role of wet (gas-rich) mergers in 
the formation of early-type galaxies, though not strongly 
constrained by the scaling relations, is limited by other 
considerations: if these systems experienced a significant 
amount of mergers in relatively recent times, most of 
these mergers must have been dry, becau se the old stel- 
lar po pulations of early- type galaxies (e.g. iThomas et al.1 
12003) are inconsistent with substantial recent star for- 
mation. It is also the case that in very massive galaxies 
most accreted cold gas is likely to be eliminated by evap- 
oration fro m the hot interstellar medium before it can 
form stars (|Nipoti fe Binnevil2007| ). 
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We emphasize that our constraints on the dry-merging 
history apply to massive early- type galaxies, of which 
the SLACS sample is representative, and are not in 
conflict with the hypothesis that exceptionally lumi- 
nous, brightest cluster galaxies (BCGs) formed via 
multiple dissipationless m ergers, the so-called galac- 
tic cannibalism scenario (jOstriker fc Tremaind 119751 : 
lHausman fc Ostrikerlll978( ). In fact, numerical simula- 
tions of dynamical-friction driven galactic cannibalism 
predict that brightest cluster galaxies must be more ex- 
tended than expected from the extrapolatio n of the size- 
luminosity correlation of smaller ga laxies (jNipoti et alj 
2003bt iRuszkowski fc Sp ringcl 2009), as observed (e.g., 



Oegerle fc Hoessellll991l;lBernar di 2009). In this respect, 



it will be very interesting to perform a systematic ex- 
ploration of the behavior of brightest cluster galaxies 
with respect to the lensing scaling relations explored 
here. Larger numbers of BCGs with lensing features and 
stellar kinematic s than those currently available (e.g., 
ISand et al.l l2008) are needed to carry out this study. 

In conclusion, we can see three main ways to make 
progress and determine whether the tightness of the lens- 
ing scaling relations is fatal to the idea of substantial 
growth of normal (i.e. non-BCG) early-type galaxies by 
dry mergers. From a theoretical point of view, fully cos- 
mological simulations are needed to correlate properties 
of the progenitors with the details of the merging process 
and verify whether the required extraordinary amount of 
fine tuning does indeed occur. Also, future studies need 



to incorporate stellar mass as a fundamental variable 
that can also be directly linked to observations. From 
an observational point of view, large and homogeneous 
samples of lenses at higher- z than SLACS are needed to 
determine how the lensing scaling relations evolve with 
redshift. 
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Seed galaxy models. 



Id 


Model 


r a /r* 


C200 




M dm /M* 


AW10 5 


7' 


A5 


A 


1.4 


7 


7.8 


10 


11 


2.02 


Al 


A 


1.4 


7 


7.8 


10 


2.2 


2.02 


B5 


B 


oo 


7 


7.8 


10 


11 


2.02 


C5 


C 


0.7 


7 


7.8 


10 


11 


2.02 


Dl 


D 


1.0 


8 


11.6 


19 


10 


2.00 



Note. — Id: name of the N-body realization. Model: name of 
the seed galaxy model. r a : anisotropy radius. C200 ; NFW con- 
centration. £ = r B /R e . M(j m : total dark-matter mass. M„: total 
stellar mass. Atot: total number of particles in the N-body model 
(stellar and dark matter particles have the same mass) . 7' : logarith- 
mic slope of the best-fit power-law total density profile (ptot oc r -7 ) 
over the radial range 0.1 < r/R c < 1. 
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TABLE 2 

Parameters of the simulations and properties of the end-products. 
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1.92 


4Alph 


2Alph 








280 


0.548 





3.95 


36.52 


0.526 


69.49 


1.95 


4Alpo 


2Alpo 





1.78 


280 


0.528 


0.130 


3.95 


35.28 


0.530 


85.17 


2.16 


4Aleh 


2Aleh 


-0.83 





280 


0.286 





3.99 


38.72 


0.613 


55.72 


1.99 


4Aleo 


2Aleo 


-0.83 


1.78 


280 


0.238 


0.120 


3.99 


38.38 


0.583 


63.14 


2.23 


4Alhh 


2Alhh 


0.83 





280 


0.678 





3.90 


32.10 


0.473 


90.19 


1.98 


4Alho 


2Alho 


0.83 


1.78 


280 


0.668 


0.122 


3.90 


29.51 


0.481 


87.52 


2.24 


4Dlph 


2Dlh 








480 


0.888 





3.89 


176.93 


0.803 


134.09 


1.70 


4Dlpo 


2Dlo 





1.78 


480 


0.864 


0.208 


3.98 


175.23 


0.804 


152.65 


1.77 


8Alph 


4Alph 








560 


0.538 





7.85 


70.78 


0.550 


160.47 


1.95 


8Alpo 


4Alpo 





1.78 


560 


0.509 


0.144 


7.82 


65.83 


0.528 


201.62 


2.18 


8Aleh 


4Aleh 


-0.83 





560 


0.066 





7.95 


76.15 


0.679 


84.12 


1.84 


8Aleo 


4Aleo 


-0.83 


1.78 


560 


0.086 


0.114 


7.95 


75.46 


0.660 


106.01 


2.16 


8Alhh 


4Alhh 


0.83 





560 


0.666 





7.63 


55.91 


0.490 


169.73 


1.96 


8Alho 


4Alho 


0.76 


0.80 


560 


0.640 


0.060 


7.60 


50.44 


0.491 


184.16 


2.13 


2A5ml0x 


A5 












1.86 


15.81 


0.587 


22.64 


2.11 


2Alml0x 


Al 












1.86 


16.09 


0.585 


23.02 


2.11 


2Alml0y 


Al 












1.86 


15.66 


0.554 


25.28 


1.84 


2Alm5 


Al 












1.76 


14.67 


0.550 


23.87 


2.03 


2DlmlO 


Dl 












1.69 


73.13 


0.896 


37.74 


1.65 


4Alml0x 


2Alml0x 












3.21 


26.99 


0.642 


32.01 


1.99 


4Alml0y 


2Alml0y 












3.61 


30.18 


0.557 


48.47 


1.94 


4Alm5 


2Alm5 












3.29 


26.69 


0.597 


38.83 


2.10 


4DlmlO 


2DlmlO 












3.07 


119.29 


0.993 


50.07 


1.66 



Note. — Id: name of the run and of the corresponding end-product. Prog.: name of the progenitor. 
E: dimensionless orbital energy. L: dimcnsionlcss orbital angular momentum modulus, do: initial 

separation. v = vo • do /do: parallel component of the initial relative velocity. Vq = ||vo X do||/do: 
orthogonal component of the initial relative velocity. M tj £ n : final stellar mass. Mdmfin'- final dark 
matter mass. <r v : virial velocity dispersion of the end-product, ryi tot'- circularized half-mass radius of 
the total density distribution of the end-product. 7': logarithmic slope of the best-fit power-law total 
density profile (ptot oc r~ 7 ) over the radial range 0.1 < r/R e < 1. Masses, lengths and velocities are 

in units of M», r* and C», respectively. The parameters E, L, do, Vq and Vq are defined only for binary 
major-merging simulations. 



